Two-Spinon and Orbital Excitations of the Spin-Peierls System TiOCl 
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We combine high-resolution resonant inelastic x-ray scattering with cluster calculations utilizing a 
recently derived effective magnetic scattering operator to analyze the polarization, excitation energy, 
and momentum-dependent excitation spectrum of the low-dimensional quantum magnet TiOCl in 
the range expected for orbital and magnetic excitations (0 — 2.5 eV). Ti 3d orbital excitations yield 
complete information on the temperature-dependent crystal-field splitting. In the spin-Peierls phase 
we observe a dispersive two-spinon excitation and estimate the inter- and intradimer magnetic 
exchange coupling from a comparison to cluster calculations. 

PACS numbers: 78.70.Ck, 71.27.+a, 71.70.Ch, 78.70.Dm 



Understanding the interplay of lattice, charge, orbital, 
and spin degrees of freedom in strongly correlated elec- 
tron systems involving frustration remains one of the 
great tasks in modern solid state physics. The role of 
electronic and magnetic excitations for, e.g., the Mott 
metal-insulator transition is a key question in this respect 
[l|. Resonant inelastic x-ray scattering (RIXS) has been 
established recently as a powerful experimental technique 
to study these phenomena in complex solids [2I413J , pro- 
moted by considerably improved energy resolution and 
detector sensitivity [1J, [T5J. While previously inelastic 
neutron scattering was the method of choice for mag- 
netic excitations, RIXS nowadays is established to give 
complementary information already from much smaller 
sample volumes (< 0.1mm 3 ). This is due to the much 
larger scattering cross section of x rays compared to neu- 
trons. Performing RIXS with photon energies close to 
the L absorption edge corresponds to processes of the 
form 2p 6 3d" ->• 2p 5 3d™ +1 -> 2p 6 3d n *. The intermediate 
state is the final state of x-ray absorption spectroscopy 
(XAS), and the RIXS final state can be an excited state 
(indicated by the asterisk) with losses due to coupling to 
all kinds of excitations. For n = 1 the initial and final 
state can be described by single-particle wave functions, 
in contrast to the intermediate state. Thus, although it is 
possible to extract the crystal-field (CF) splitting directly 
from orbital (dd) excitations, a many-body approach is 
necessary to simulate RIXS spectra. 

In this Letter, we analyze the temperature, excita- 
tion energy, polarization, and momentum dependence 
of the magnetic and orbital excitations of the quasi- 
onc-dimcnsional spin-Peierls system TiOCl measured by 
high- resolution RIXS. The data allow us to extract pre- 
cise values for the CF splitting in the room temperature 
(RT) as well as the spin-Peierls phase. The polarization 
dependence of the CF levels is in very good agreement 
with cluster calculations based on density-functional the- 



ory within the local-density approximation. From com- 
parison to cluster calculations for a dimcrized spin chain 
we identify a low-energy feature as a two-spinon exci- 
tation, amplifying the quasi-one-dimcnsional nature of 
TiOCl. The predictive power of a recently derived effec- 
tive scattering operator for magnetic excitations in RIXS 
is evidenced for the first time by simulating the intensity 
of this excitation for several different momenta in relation 
to the d-d excitation strength. Furthermore, by includ- 
ing the dimerization into our theory such that the spin 
gap found in other experiments [16[ is reproduced we can 
infer the ratio of the inter- and intradimer magnetic ex- 
change coupling. 

TiOCl is a layered system of buckled Ti-0 bilayers 
sandwiched by CI ions. Such layers are stacked along 
the c axis with only weak van der Waals forces present in 
between, allowing for easy cleavage. The structure thus 
effectively confines electronic and magnetic interactions 
to the ab plane. In this plane, Ti ions form a staggered 
triangular lattice which introduces magnetic frustration. 
Locally, each Ti ion is coordinated by four O and two 
CI ions which form a strongly distorted octahedron with 
a local C'2v point-group symmetry in the orthorhombic 
crystal structure. This results in a completely lifted de- 
generacy of the 3d orbitals, i.e., five energetically well 
separated CF levels, and thus a quenching of the orbital 
degrees of freedom 17j • The Ti ions are in a 3d 1 configu- 
ration and the overlap of occupied d x 2_ y 2 orbitals along b 
provides a one-dimensional (ID) hopping path, leading to 
(quasi-) ID behavior in several aspects [lj|. For example, 
below T cl = 67 K the system shows a dimerization of Ti 
ions along b into S = singlets, i.e., a spin-Peierls phase, 
with a huge spin gap of approximately 430 K = 37meV 
[161 ]. At higher temperatures the magnetic susceptibil- 
ity corresponds to that of a ID Heisenbcrg chain, with a 
nearest-neighbor exchange coupling J « 660 K= 57meV 
[191 ] . The spin-Peierls ground-state below 67 K is of con- 
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FIG. 1: (Color online) RIXS spectra at T = 300 K (light red) 
and T = 36 K (dark blue) (hv = 454.2 eV, a polarization, 
ac scattering; see text). Black solid curves are Gaussian fits 
to the low-T data. Also, the orbital character of the CF ex- 
citations is indicated. Insets: Experimental RIXS geometry. 
Left: (j> = 90° geometry to define a (out-of-plane) and n (in- 
plane) polarizations. Right: Definition of q (brown arrow) 
and Q|| (red arrow), q was fixed by the constant 4> — 50° in 
our experiments, and qu was varied by rotating the sample 
around the polar axis, i.e., changing 9. 



siderable interest since it is reached via an incommensu- 
rately modulated intermediate phase below T C 2 = 91 K 
[20l ] . Frustrated interchain interactions are widely ac- 
cepted as the driving force for this unconventional be- 
havior (cf. Ref. [21| and references therein). 

Samples were grown by a chemical vapor transport 
method [22j, with a typical size of (1 x 2 x 0.1) mm 3 . 
Measurements were conducted at the SAXES endstation 
UGH of the ADRESS beamline (spot size 52x4^m 2 ) at 
the Swiss Light Source (Paul Scherrer Institut, Switzer- 
land) using soft x-ray RIXS with an energy resolution 
of 80meV or better, as determined from the full width 
at half maximum of the elastic line of a carbon powder 
sample at the relevant excitation energies. As shown in 
the insets of Fig. Q] the endstation allows different angles 
4> between the incoming and outgoing beam, as well as 
a rotation of the sample around the polar axis for the 
momentum-dependent measurements. For = 90° and 
7r-polarized (i.e., within the scattering plane) incoming 
light the polarization of the outgoing beam must nec- 
essarily be rotated. For er-polarized (i.e., perpendicu- 
lar to the scattering plane) incoming light this need not 
be the case. This leads to well-defined polarization ef- 
fects even without detecting the polarization of the out- 
going photons. We used <j> = 90° and a grazing in- 
cidence angle 9 = 20° for all RIXS measurements ex- 
cept those to investigate the momentum dependence. In 
the latter case, only the projection of the momentum 
transfer q = k m — k oxlt on the sample surface (q\\, see 
right inset of Fig. [1]) is of importance due to the quasi- 
1D nature of TiOCl. Because of the bulk sensitivity of 



RIXS governed by the attenuation length of the x rays 
(> 100 A), it was sufficient to cleave samples in air us- 
ing Scotch tape and then transfer them into the vacuum 
(< 10 _9 mbar). Acquisition times per spectrum were in 
the range 90 — 200 min. To simulate the dd excitations 
we used a (paramagnetic) TiGjC^ cluster within multi- 
plet ligand-field theory, a powerful and well-established 
many-body approach to x-ray absorption and emission 
processes, especially in strongly correlated electron sys- 
tems |23j . Hopping parameters and Slater integrals were 
obtained from local-density approximation calculations 
followed by a downfolding procedure [24|, |25J ■ The clus- 
ter calculations account for the CF splitting as well as 
the covalence between the Ti 3d and O 2p orbitals. The 
magnetic excitations have been calculated from two sepa- 
rate clusters utilizing a recently derived effective scatter- 
ing operator for RIXS of collective magnetic excitations 
[101 ] . The x-ray absorption spectra that enter within this 
theory as the energy-dependent resonant enhancement 
prefactors were obtained from the same cluster as the 
dd excitations. For the two-spinon excitations we used a 
cluster of 28 sites for a ID dimerized Hciscnbcrg chain 
with nearest-neighbor exchange. 

Figure [T] shows typical RIXS spectra at T = 300 K (red 
curve) and at T — 36 K (blue curve), i.e., in the spin- 
Pcicrls phase, taken at a photon energy hv = 454.2 eV 
in a polarization with ac being the scattering plane ( "ac 
scattering"). In both cases, five peaks are observable, 
corresponding to the on-site dd excitations. For the first 
time, all CF excitations of TiOCl are mapped directly in 
one experiment, as in other (optical) probes some of them 
are forbidden due to dipole selection rules [26|. Besides 
the elastic line at zero ener gy loss four peaks appear at 
0.34, 0.62, 1.48, and 2.08 eV [3lj|, as determined from fit- 
ting Gaussian peaks to the low-T data. At 300 K, they 
lie at 0.36, 0.59, 1.45, and 2.05 eV, respectively, due to T- 
dependent changes in the distortion of the Ti04Cl2 oc- 
tahedra predicted by cluster calculations [26[ . 

Figure [5{a) shows RIXS spectra (solid red curves) at 
different photon energies across the Ti L edge (T = 
300 K). Also shown are results from cluster calculations 
(dashed blue curves), including 60meV resolution and 
200 meV (final-state) lifetime broadening as inferred from 
experiment. The elastic line has zero weight in the sim- 
ulations (except for Bragg reflections). This is expected 
for the perfectly ordered, infinite crystal assumed in our 
calculations. In experiment, surface reflectivity and scat- 
tering off crystal imperfections always produce a (quasi- 
elastic peak. The overall agreement between theory and 
experiment is very good, and aside from sample quality 
the remaining discrepancies might be connected to un- 
certainties concerning the exact excitation energy. This 
is illustrated by Fig. EJb) which shows absorption spectra 
acquired in total electron yield, with the (main) polar- 
ization [321 ] and the used photon energies indicated. The 
XAS consists of a complicated multiplet structure, and 
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FIG. 2: (Color online) (a) Excitation-energy dependence of RIXS spectra at T — 300 K in a polarization and ac scattering. 
Corresponding excitation energies are indicated, (b) Ti L edge XAS spectra with main polarizations along the different crystal 
axes both measured (solid curves) and calculated (dashed curves) for the geometry of panel (a) . Photon energies from (a) are 
marked by vertical lines. (c),(d) a- and 7r-polarized RIXS spectra for different photon energies in ac (c) and be (d) scattering 
from experiments (solid curves) and multiplet calculations (broken curves). 



being off in photon energy by a few meV can lead to 
largely different intensities for the involved dipole tran- 
sitions. Since XAS is the intermediate state of RIXS 
the agreement between simulated and measured RIXS 
spectra cannot be better than for the absorption spec- 
tra. Note that the intermediate state broadening for 
the RIXS spectra is equivalent to the broadening used 
for the calculation of the XAS spectra. The core-hole 
broadening varies between 200 meV for the low-energy 
states (450 eV) to 600 meV for the high energy states 
(465 eV). The same also holds for the spectra shown in 
Figs. Wis) and (d), where the effects of different polar- 
ization and scattering geometries are compared at two 
excitation energies. Again, although some discrepancies 
remain there is notable agreement between theory and 
experiment. Note that in the experimental data the elas- 
tic line is strongly suppressed in it polarization 27], |28| . 



A feature barely visible in the experimental spectra 
shown up to now is a low-energy excitation (LEE) cen- 
tered at roughly 0.1 eV, which is not captured by the 
paramagnetic cluster calculations. Figure G3a) shows 
momentum-dependent RIXS spectra within an energy 
window containing only the two lowest dd excitations 
and the elastic line measured at T — 15 K for different 
q\ | along b (it polarization, be scattering) . Since for some 
samples the LEE feature started to disappear after 5 min 
of beam exposure, the sample spot was changed every 
5 min by 10 fan. Note, however, that also for such sam- 
ples the behavior of the LEE completely agrees with the 
data shown here, taken on samples, which did not show 
any degradation of the spectra up to 20 min of irradia- 
tion. Simulated spectra in Fig.O^a) consist of a superpo- 
sition of three parts. In addition to dd excitations with 
a reduced lifetime broadening of 150 meV due to the low 
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FIG. 3: (Color online) (a) RIXS spectra taken at hu = 
456.8 eV in n polarization and be scattering (T = 15 K), for 
q\\ corresponding to (from top to bottom) 36%, 31%, 25%, 
19%, 11%, 0%, and -11% of the FY path in the BZ. The LEE 
is highlighted by the shaded area, (b) Two-spinon dynamical 
spin structure factor for a 28 site dimerized Heisenberg model. 
The dispersion of the LEE is indicated, obtained from fitting 
the data in panel (a) (blue square markers). At mi =0 the 
spin gap extracted from NMR [1(| is indicated. Inset: LEE 
spectral weight (i.e., experimental peak areas). 



measuring temperature, spinon calculations for the LEE 
based on the above-mentioned effective operator [10| , and 
a Gaussian to mimic the elastic peak are incorporated. 
Since at hv = 460 eV one has q « 0.413 A -1 and the Bril- 
louin zone (BZ) boundary along the ID direction (the Y 
point) lies at 0.94A -1 wc could cover approximately 1/3 



of the path TY with q\ j . 

The LEE is most pronounced in the spectrum with 
the highest q l{ = 0.338A" 1 [top of Fig. EJa)]. Towards 
the r point it rapidly loses intensity, and a dispersion 
from ~ 110 meV to ~ 80meV is observed. Its nature 
becomes evident from looking at the tw o-sp inon dynam- 
ical structure factor for a ID chain [29|, |30|. It describes 
a continuum of excited states with two sinelike (quasi- 
)boundaries wz,(g) and uju(q), and its spectral weight 
Ilj U () ^(9> a - , )dw goes linearly to zero towards the zone 
center. To take dimerization into account we calculated 
the dynamical structure factor for a 28 site Heisenberg 
spin chain with intra- and interdimer nearest-neighbor 
interactions J and J', respectively. As can be seen in 
Fig-E^b) a spin gap opens at the zone center and the zone 
boundary. The intensity, however, still goes linearly to 
zero towards the zone center as for the undimcrized chain. 
Note that in the experiment the in-plane momentum and 
the polarization are changed simultaneously between dif- 
ferent spectra in Fig.|3|a) due to the experimental setup. 
This has been accounted for in the calculations, but the 
dependence of the spectral weight on the polarization for 
this specific case is only of the order of 10%. Conse- 
quently, there is good agreement for the intensity ratios 
of both the spin and the orbital excitations. The experi- 
mental dispersion (blue square markers) and the spin gap 
from NMR (red dot, 37meV) [16[ are best reproduced 
for exchange couplings of J = 90mcV and J' = 0.75 J, 
respectively, i.e., their mean value is in reasonable agree- 
ment with the value of 57 meV for the undimcrized chain. 
The inset shows that the experimental LEE weight goes 
linearly to zero, evidencing its two-spinon character. An 
analogous behavior of magnetic features in RIXS has pre- 
viously been observed e.g. on La2CuC>4 
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is the first noncuprate in which such an excitation could 
be identified. 

In principle, dd excitations might also be dispersive. 
The orbital exchange energy is of the same order as the 
magnetic exchange and smaller than their lifetime broad- 
ening. Since in RIXS the probability for flipping an or- 
bital or a spin is roughly the same the experimental exci- 
tations likely involve coupled orbital and magnon excita- 
tions, eventually with additional damping from phonons. 
Looking back at Fig. |3ja) one notices that at large g||, 
e.g., the d yz excitation is significantly broader than the 
local cluster calculation. A detailed analysis of this obser- 
vation, however, is out of the scope of the present study. 

In conclusion, by overcoming the selection-rule limita- 
tions of optical spectroscopy using RIXS we have directly 
obtained reliable values for the CF splitting from the or- 
bital (dd) excitations of TiOCl. Their temperature, ex- 
citation energy, and polarization dependence is in agree- 
ment with results from cluster calculations. Evaluating 
the momentum-dependent behavior of a low-energy peak 



observed in the spin-Peierls phase on the energy scale of 
the magnetic exchange coupling J, we can identify this 
peak as a two-spinon excitation. Using a new effective 
operator for magnetic excitations in RIXS in combination 
with finite size cluster calculations its dispersion and in- 
tensity is accounted for. Moreover, we find robust values 
for the inter- and intra-dimer exchange coupling. 
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